The Limits of Extended Quintessence 
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We use a low redshift expansion of the cosmological equations of extended (scalar-tensor) 
quintessence to divide the observable Hubble history parameter space in four sectors: A forbid- 
den sector I where the scalar field of the theory becomes imaginary (the kinetic term becomes 
negative), a forbidden sector II where the scalar field rolls up (instead of down) its potential, an 
allowed 'freezing' quintessence sector III where the scalar field is currently decelerating down its 
potential towards freezing and an allowed 'thawing' sector IV where the scalar field is currently 
accelerating down its potential. The dividing lines between the sectors depend sensitively on the 
time derivatives of the Newton's constant G over powers of the Hubble parameter. For minimally 
coupled quintessence which appears as a special case for a constant G our results are consistent 
with previous studies. Observable parameter contours based on current data (SNLS dataset) 
are also constructed on top of the sectors, for a prior of Qm = 0.24. By demanding that the ob- 
served 2a parameter contours do not lie entirely in the forbidden sectors we derive stringent 
constraints on the current second time derivative of Newton's constant G. In particular we find 
^ > —1.91 Hq — —2 X I0~^''h^ yrs~^ at the 2g level which is complementary to solar system tests 
which constrain only the first derivative of G as |^1 < lQ~^^yrs~^ at la. 

PACS numbers; 98.80.Es,98.65.Dx,98.62.Sb 



I. INTRODUCTION 

Cosmological observations based mainlyon type la su- 
pernovae (Snia) standard candles [J and also on 

standard rulers 0, [1] have provided a fairly accurate 
form of the Hubble parameter H{z) = ^ as a function 
of redshift z in the redshift range < z < 1.7 ^ and 
beyond [6]. This form indicates that despite the attrac- 
tive gravitational properties of matter, the universe has 
entered a phase of accelerated expansion at a redshift 
z ~ 0.5. It is therefore clear that the simplest cosmolog- 
ical model, where the universe is dominated by matter 
and its dynamics is determined by general relativity is 
ruled out at several a The central current questions 
in cosmology research are the following: 

• What theoretical models are consistent with the 
currently detected form of H{z)l 

• What are the generic predictions of these models 
with respect to H{z) so that they can be ruled 
out or confirmed by more detailed observations of 

H{z)l 

In a class of approaches the required gravitational 
properties of dark energy (see 0, H, H, [io| for recent re- 
views) needed to induce the accelerating expansion are 
well described by its equation of state w{z) = . In 
this case the simplest model consistent with the currently 
detected form of H{z) is the flat KCDM model. Accord- 
ing to this model, the universe is flat and its evolution 
is determined by general relativity with a cosmological 
constant through the Friedman equation 



where fion 



is the current matter density normal- 



H{zf = Hl{no^{i + zf + 



(1.1) 



ized on the critical density for flatness pc and Oa = 
1 — r^om is a constant density due to the cosmological 
constant. The main advantages of this model are simplic- 
ity and predictability: It has a single free parameter and 
it can be definitively ruled out by future observations. 
Its disadvantages are lack of theoretical motivation and 
fine tuning: There is no physically motivated theoretical 
model predicting generically a cosmological constant at 
the observed value. This value is 120 orders of magnitude 
smaller than its theoretically expected value [111, ■ 

Attempts to replace the cosmological constant by a 
minimally coupled dynamical scalar field (minimally cou- 
pled quintessence (MCQ) [11, Hi, Ull ) have led to mod- 
els with a vastly larger number of parameters fueled by 
the arbitrariness of the scalar field potential. Despite 
this arbitrariness and vast parameter space, quintessence 
models are generically constrained to predicting a lim- 
ited range of functional forms for H{z) [H, [l3|- This is 
a welcome feature which provides ways to either rule out 
or confirm this class of theories. 

The allowed functional space oi H{z) can be further in- 
creased by considering models based on extensions of gen- 
eral relativity such as braneworlds ITol. [20l [2ll . f(R) 
theories [l^ or scalar-tensor theories [l6l nam 

(ex- 
tended quintessence (EXQ) [2J]) where the accelerated 
expansion of the universe is provided by a non-minimally 
coupled scalar field. This class of theories is strongly mo- 
tivated theoretically as it is predicted by all theories that 
attempt to quantize gravity and unify it with the other 
interactions. On the other hand, its parameter space is 
even larger than the corresponding space of MCQ since 
the later is a special case of EXQ. Local (eg solar system) 
gravitational experiments and cosmological observations 
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constrain the allowed parameter space to be close to gen- 
eral relativity. Despite of these constraints however, the 
allowed by EXQ functional forms oiH{z) are significantly 
more than those allowed by MCQ. The detailed identi- 
fication of the forbidden H{z) functional forms for both 
MCQ and EXQ is particularly important since it may 
allow future observations determining H{z) to rule out 
one or both of these theories. 

Previous studies [1^ have mainly focused on the H{z) 
limits of MCQ using a combination of plausibility ar- 
guments and numerical simulations of several classes of 
potentials. The low redshift H{z) parameter space was 
divided in three sectors: a forbidden sector which could 
not correspond to any plausible quintessence model, a 
sector corresponding to the freezing quintessence scenario 
and a sector corresponding to the thawing quintessence 
scenario. In the freezing quintessence models, a field $ 
which was already rolling towards its potential minimum 
prior to the onset of acceleration slows down ($ < 0) 
and creeps to a halt (freezes) mimicking a cosmological 
constant as it comes to dominate the universe. In the 
thawing quintessence models, the field has been initially 
halted by Hubble damping at a value displaced from its 
minimum until recently when it 'thaws' and starts to roll 
down to the minimum ($ > 0). 

Here we extend these studies to the case of EXQ. In- 
stead of using numerical simulations however, applied 
to specific potential classes, we use generic arguments 
demanding only the internal consistency of the theory. 
Thus our 'forbidden' sector when reduced to MCQ is 
smaller but more generic (applicable to a more general 
class of models) than that of Ref. [25| (see also Refs. 

The size and location of the sectors of the low z, H{z) 
parameter space, depends sensitively on the assumed 
current time derivatives of the Newton's constant G{t) 
and reduce to well known results in the MCQ limit of 
G{t) = Gq = const. Therefore an interesting interplay 
develops between local gravitational experimental con- 
straints of the current time derivatives of G{t) (eg ^ or 
^) and cosmological observations of H{z) at low red- 
shifts. For example, a constraint on ^ — ^ from local 
gravitational experiments defines the forbidden sector in 
the low z expansion coefficients of H{zY in the context 
of EXQ. If such coefficients are measured to be in the 
forbidden sector by cosmological observations then EXQ 
could be ruled out. Alternatively, if such coefficients are 
measured to be in the forbidden sector for MCQ but in 
the allowed sector of EXQ (either 'freezing' or 'thaw- 
ing') then this would rule out MCQ in favor of EXQ. As 
shown in what follows, current observational constraints 
on H{z) imply significant overlap with the allowed sectors 
of both MCQ and EXQ. This however may well change 
in the near future with more accurate determinations of 
H{z) and the time derivatives of Gq. 



II. THE BOUNDARIES OF EXTENDED 
QUINTESSENCE 

Extended quintessence is based on the simplest but 
very general (given its simplicity) extension of general 
relativity: Scalar-Tensor theories. In these theories New- 
ton's constant obtains dynamical properties expressed 
through the potential F{^). The dynamics are deter- 
mined by the Lagrangian density [H, 111] 



, (2.1) 

where Cm[ipm', g^iv] represents matter fields approximated 
by a pressureless perfect fluid. The function F{^) is ob- 
servationally constrained as follows: 

• F{^) > so that gravitons carry positive 
energy [23j. 

• -pr (^)^ |2=o < 3:10"^ from solar system observa- 
tions [28]. 

In such a model the effective Newton's constant for the 
attraction between two test masses is given by 

\ 2 



Geffit) 



1 F{t) + 2{%{t)y 1 



F{t) 



G{t) (2.2) 



where the approximation of equation (j2.2|) is valid at low 
redshifts. Assuming a homogeneous and varying the 
action corresponding to (|2.ip in a background of a flat 
FRW metric, we flnd the coupled system of generalized 
Friedman equations 



iFH^ 
-2FH 
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-<i>^ - 3HF 
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(2.3) 
(2.4) 



where we have assumed the presence of a perfect fluid 
{p = Pm,P — 0) playing the role of matter fields. Ex- 
pressing in terms of redshift and eliminating the potential 
U from equations (USD, (EH) we find 



-F" 



{InH)' 
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2^^F 



-3{l + z)fl, 
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Fo 



(2.5) 



where ' denotes derivative with respect to redshift and Fq 
is set to 1 in units of ■ 



SttGc 



and corresponds to the present 
value of F. Alternatively, expressing in terms of redshift 
and eliminating the kinetic term from equations 
we find 



U 



(l + z)^g^ 
2 
6 

(l-t-z)2 " 



(IniJ)' 
(IniJ)' 



F' 



- 3(1 + .) (f 



l + z 
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(2.6) 



3 



We now wish to explore the observational consequences 
that emerge from the following generic inequalities antic- 
ipated on a purely theoretical level 



> 
U'{z) > 



(freezing) 
(thawing) 



(2.7) 
(2.8) 
(2.9) 
(2.10) 



The inequality (|2.7p is generic and merely states that 
the scalar field in scalar-tensor theories is real as it is 
directly connected to an observable quantity (Newton's 
constant). The inequality (|2.8p is also very general as it 
merely states that the scalar field rolls down (not up) its 
potential. This inequality is not as generic as (|2.7p since 
it implicitly assumes a monotonic potential. Finally, 
the inequality (|2.9p ( (|2.10p ) denotes a scalar field which 
decelerates (accelerates) as it rolls down its potential 
thus corresponding to a freezing (thawing) quintessence 
model. 

Since we are interested in the observational conse- 
quences of equations (|2.7p - (|2.10p at low redshifts, we 
consider expansions of equations (|2.5p and (|2.6p around 
z = expanding F{z), H{z)^ U(z), $(z) as follows 



F(z) 
H(z)' 
U(z) 
$(z) 



= l + Fiz + F2Z^ + . 
= l + hiz + h2Z^ + . 

= l + UiZ + U2Z^ + 



(2.11) 
(2.12) 
(2.13) 
(2.14) 



where we have implicitly normalized over Fq , Hq , Uq and 

$0- 

It is straightforward to connect the expansion coeffi- 
cients Fi of equation (j2.1ip with the current time deriva- 
tives of G(t) using equation (|2.2p and the time-redshift 
relation 



dt 

dz 



1 



H{z)(l + z) 
For example for Fi we have 



' F^ dz Goi/o 
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(2.15) 



(2.16) 



where the subscript o denotes the present time and Hq 
10""'^'' h yrs~^. Similarly, for F2 we find 



F2 = gi(gi 
where we have defined 

9n = 



hi 



92 

2 



GqHq 



(2.17) 



(2.18) 



with the superscript ^"^ denoting the n*'* time derivative. 
We may now substitute the expansions (|2.1ip . (|2.12p . 
(I2.14P in equation (|2.5p replacing the coefficients Fi by 



the appropriate combination of Equating terms order 
by order in z and ignoring terms proportional to gi due 
to solar system constraints which imply [sO, [HI 



|5i| < lO-'hjrs-^H-' 



IQ-^h-' < 1 



(2.19) 



we find for the zeroth and first order in 



hi - Sf^om + .g2 = > (2.20) 
-hi(l + hi) + 2h2 - 3f^om(l - hi)- 
-52(1 + hi) - 33 - 4$i$2 = (*")'(^ = 0) (2.21) 



The inequality (|2.20p defines a forbidden sector for ex- 
tended quintessence for each value of 32- For (?2 = 
this reduces to the well known result that MCQ can not 
cross the phantom divide line (see equation (|2.40p be- 
low). Equation f^?^ can be used along with (PTU|) 
to divide the allowed (hi,h2) parameter sector into a 
freezing quintessence sector ($'(z)^)' > and a thawing 
quintessence sector ($'(2;)^)' < for each set of (52,33)- 

Unfortunately, solar system gravitational experiments 
have so far provided constraints for gi [s^] (equation 
(|2.19p ) but not for gi with i > 2. This lack of con- 
straints is not due to lack of observational data quality 
but simply due to the fact that existing codes have pa- 
rameterized G(t) in the simplest possible way ie as a lin- 
ear function of t. It is therefore straightforward to extend 
this parameterizations to include more parameters thus 
obtaining constraints of Qj with i > 2. Such an analysis 
is currently in progress [32l | but even before the results 
become available we can use some heuristic arguments 
to estimate the order of magnitude of the expected con- 
straints on gi (i > 2) given the current solar system data. 

Current solar system gra vity experiments are utiliz- 
ing lunar laser ranging [Slj and high precision planet 
ephemerides data [30| to compare the trajectories of ce- 
lestial bodies with those predicted by general relativity. 
The possible deviation from the general relativity pre- 
dictions is parameterized [31] using the PPN parameters 
f3 and 7, the first current derivative of Newton's con- 
stant ^ and a Yukawa coupling correction to Newton's 
inverse square law. These experiments have been col- 
lecting data for a time At of several decades (sol. l3lj | ie 
At = O(lOOyrs). The current la constraint [3^ of 



1^1 ^ 1-0.2 ±0.5| xlO^^yrs-i < 10"" 
Go 



yrs 



(2.22) 



implies an upper bound on the total variation over 
the time-scale At of approximately 

,AG, ,Go 



Go Go 



\At < 10' 



(2.23) 



The same bound is obtained by considering the relative 
error in the orbital periods T of the Earth and other 
planets which at 1 cr is [s^l 



AT 



< 10 



-12 



(2.24) 
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FIG. 1: The hi — /12 sectors of EXQ: Sector I is the forbidden sector where the scalar field becomes imaginary. Sector II is 
also forbidden and corresponds to a scalar field that rolls up its potential. Sector III corresponds to 'freezing' EXQ where 
the field is decelerating down its potential towards freezing. Sector IV corresponds to 'thawing' EXQ where the scalar field is 
accelerating down its potential. The forbidden sector I shrinks for a Newton's constant G that decreases with time (172 > 0) 
while the 'freezing sector' expands (Fig. la). The boundaries of the sectors are provided by (|2.2U|l . (|2.2ip and (|2.31|l . The la, 
2(7 contours obtained from the SNLS SnIa dataset for flom = 0.24 using the GPL parameterization, are also shown. 



Given that the Keplerian orbital period is 



we find 



AT 
I — I 



Gra 



1 AG 

2' Go ' 



< 10" 



(2.25) 



(2.26) 



in rough agreement with (|2.23p . 

Using the upper bound p.23p and attributing any vari- 
ation of G to a term quadratic in At we get 



l^l-l^l(At)^<10-" 



which imphes that 

1^1 < lo-^V^-^ 



1.92 1 < lo^/i- 



(2.27) 



(2.28) 



giving a rough order of magnitude estimate of the up- 
coming constraints on 32- Preliminary results from the 
analysis of solar system data indicate that [111 



Go 

Go 



(4±5) • lO-^^yrs-2 



(2.29) 



which is not far off the rough estimate of equation (|2.28p . 
Generalizing the above arguments to arbitrary order in 
A< we find 



\9n\ < 10 



8n — 11 u—n 



(2.30) 



which implies that current solar system tests can not con- 
strain gn in any cosmologically useful way for n > 2. 



We now return to the constraint equations (|2.20p and 
(I2.2ip and supplement them by the constraint obtained 
from the inequality (|2.8[) . Using the expansions (|2.1ip . 
(|2.12p and (|2.13p in equation (|2.6p and equating terms of 
first order in z we find 



Ui = U'{z ^0)^ - {5hi 



2ho 



gr^om + 5^2 + 93) > 

(2.31) 

where as usual we have ignored terms proportional to gi 
due to ([2A9| . We now may use jMQl), ([rgTjl and ([23T|) 
to define the following sectors in the hi — /i2 parameter 
space for fixed gi, g2- 



• A forbidden sector I where the inequality (|2.20p is 
violated. 



• A forbidden sector II where the inequality (I2.3ip is 
violated but not the inequality (|2.20p . 

• An allowed sector III of freezing EXQ where the in- 
equalities (|2.7p . (|2.8p and (|2.9p are respected while 
((TTU)) is violated. 

• An allowed sector IV of thawing EXQ where the in- 
equalities (|2.7p , (|2.8p and (|2.10p are respected while 
HiM is violated. 



These sectors are shown in Fig. 1 for Jlom = 0.24[6[, 
gs = and three values of g2'. g2 — 1.97, ^2 = and 
92 = -1.91. 

It would be interesting to superpose on these sectors of 
Fig. 1 the parameter contours obtained by fitting to 
the SnIa data. It is not legitimate however to extrapolate 
the expansions in powers of z out to z = 0(1) where the 
SnIa data extend. We thus use an extrapolationjs^, [s^ 
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FIG. 2: The sectors I-IV of Fig. 1 mapped on the wg — wi parameter space. The aUowed sectors extend to uiq < — 1 for g2 > 0. 



with several useful properties first proposed by Cheval- 
lier and Polarskijs^l and later investigated in more detail 
by LinderfsH. We superpose on the sectors of Fig. 1, the 
la and 2a hi — /i2 contours obtained by fitting the 
Chevalier-Polarski-Linder (CPL) [s^, [sH] parametriza- 
tion 

+ (1 - 17o™)(l + z)='(i+'""+"'i)eTS^f ](2.32) 
to the SNLS dataset [l] with the same ilom 

prior (r^om — 

0.24). Using equation (|2.32p along with the expansion 
(|2.12p it is easy to show that 

hi = 3(1 + wo- floniWo) (2.33) 
3 

/l2 = -(2 + 5wo(l-f^Om) + 

+ {I - nQ^)i3wl + wi)) (2.34) 

so that the standard contours in the wq — wi space 
(see eg [36[) can be easily translated to the — /12 space 
of Fig. 1. The CPL parametrization (|2.32p is constructed 
so that the parameter w{z) defined as 



to be of the form 



(l + z)^ 



dlnH 



1 



iLa. 



n^^ii + zf 



w{z) — Wq+Wi 



1 + 2 



(2.35) 



(2.36) 



which has useful and physically motivated properties [35|. 
The parameter w{z) is particularly useful and physically 
relevant in the MCQ limit of gi 0. In that limit w{z) 
becomes the MCQ equation of state parameter ie 



w{z) 



PMCQ _ - 
PMCQ + f7($) 



(2.37) 



as may be shown from the MCQ Friedman equations. 
The values of §2 used in Figs la and Ic were motivated 
by demanding minimum and maximum overlap of the 
forbidden sector I with the 2a hi — /12 contour. 

The following comments can be made with respect to 
Fig. 1: 

• For values oigi < —1.91, the 2a parameter contour 
obtained from SNLS lie entirely in the forbidden 
sector. This bound is independent of 93 which does 
not enter in the inequality (j2.20p . Thus in the con- 
text of EXQ the constraint on g2 obtained by the 
SNLS data at 2a level is 



92 



GqHq 



> -1.91 



(2.38) 



Notice the dramatic improvement of this constraint 
(with respect to the lower bound) compared to the 
anticipated constraint of —10^ < <?2 < 10^ antici- 
pated from solar system data (equation (|2.28p )! 

• The parameter sector III of freezing quintessence 
is significantly smaller than sector IV of thawing 
quintessence and the difference is more prominent 
for smaller g2. 

• For (72 > the allowed parameter space increases 
significantly compared to MCQ (32 = 0). There- 
fore, if future cosmological observations show pref- 
erence to the forbidden sectors I or II of Fig. lb 
(MCQ) this could be interpreted as evidence for 
EXQ with 52 > (see also comments in Refs. 
0,^,133) 

Even though the plots of Fig. 1 capture the full physi- 
cal content of our results it is useful to express the sectors 
I-IV in terms of parameter pairs other than hi — /i2 which 
are more common in the literature. Such parameters are 
the expansion coefficients Wi of w{z) 



w{z) = Wq + WiZ + W2z'^ + 



(2.39) 
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which is connected to H{z) via equation (|2.35[) . By ex- 
panding both sides of equation (|2.35p using equations 
()2.12p and (|2.39p we may express hi in terms of Wi thus 
rederiving equations (|2.33p and (|2.34p for i — and i — 1. 
The resuh is identical since the wq, wi expansion coeffi- 
cients of equation (I2.39P coincide with the wq, wi coef- 
ficients of the CPL parametrization for a w{z) given by 
equation (j2.36p . The advantage of using the parameters 
Wq — wi instead of hi — /12 is that they provide better 
contact with previous studies and can illustrate clearly 
the fact that a (72 > can provide a phantom behavior 
Wq < —1 and crossing of the phantom divide line w = — 1 
in EXQ models. Using equations (|2.33p and (|2.34p we 
may express the constraint equations ()2.20p . ()2.2ip and 
()2.3ip that define the sectors I-IV in Fig. 1 in terms of 
Wq — Wi. The resulting equations are 

3{l ~ nom){l + wo) + 92 = > (2.40) 

3(1 - f2om)((l + Wo)(3iloinWo - 2) + Wi)- 

-(72(4 + 3(1 - nom)wo) ' .93 = (*'')'(2 = 0) (2.41) 
1(1 - 17o,„)(3(l - wl) - u;i) + §52 + f > (2.42) 

In the MCQ limit (32 ^ 0, 53 ^ 0) equation ([2:421) 
has also been obtained in Ref. as a generic limit of 
MCQ. Using now ((2:401) - ((2:421) along with the wn_- wi 
contours obtained from the SNLS dataset P, [3^ we 
construct Fig. 2 which is a mapping of Fig. 1 on the 
Wo — Wl parameter space. An interesting point of Fig. 2 
is that for 52 > (Fig. 2a) the forbidden sector I shrinks 
significantly compared to MCQ (Fig. 2b) and allows for 
a, Wo < —1. 

A final set of parameters we consider is the set of the 
expansion coefficients of the luminosity distance d^^z) 
which in a flat universe is connected to the Hubble ex- 
pansion history H{z) as 

Expanding dL{z) as 

dL{z)^ z + dL2Z^ + dL3.z^ + ... (2.44) 

and using the expansion (I2.12p in equation (|2.43p we may 
express the coefficients hi — /12 in terms of ^^2 — c?L3 as 
follows 

hi = 4(1 -di2) (2.45) 
/i2 = 6(1 - 3rfL2 + 2^2 - ^^^3) (2.46) 

Substituting now equations (|2.45p . (|2.46p in equations 
()2.20p . (|2.2ip and ()2.3ip we obtain the sector equations 
in dL2 — dL3 space as 

4(1 - dL2) ~ 3f^0m + .92 = > (2.47) 
4di2(2dL2 +92- 3f^0m) + Q^On,- 

-8 - 12di3 - 552 - .93 - (3''')'(2 - 0) (2.48) 

4 + 8dL2 - 12^1,2 + 6dL3 - f r!om + 

+ i<72 + f > (2.49) 



Using now equations ((2.45p . ((2.46P to translate the hi — h2 
contours to the dL2 ~ c?L3 parameter space and equa- 
tions (|2.47p - (|2.49p to construct the sectors I-IV in the 
dL2 — dL3 space we obtain Fig. 3. The advantage of Fig. 
3 compared to Figs. 1 and 2 is that it refers to the pa- 
rameters dL2 — dL3 which are directly observable through 
the luminosity distances of Snia without the need of any 
differentiation. 



III. CONCLUSION-DISCUSSION-OUTLOOK 

We have used generic theoretically motivated inequal- 
ities to investigate the space of observable cosmological 
expansion parameters that admits viable MCQ and EXQ 
theoretical models. Our inequalities are generic in the 
sense that they are independent of the specific features of 
any scalar field potential (eg scale, tracking behavior etc) 
and they only require that the models are internally con- 
sistent. The derived forbidden sectors which violate the 
above inequalities already have significant overlap with 
the parameter space which is consistent with observa- 
tions at the 2(7 level. This overlap which depends on the 
time derivatives of the Newton's constant G has lead to a 
useful constraint to the second derivative of G (equation 
(|2.38p ) which is significantly more stringent compared to 
the corresponding constraint ((2.28P anticipated from so- 
lar system gravity experiments . 

An important reason that limits the observable param- 
eter space consistent with MCQ and EXQ is the fact that 
the scalar field potential energy can induce accelerating 
expansion but not beyond the limit corresponding to the 
cosmological constant {w{z) = —1) obtained when the 
field's evolution is frozen. Additional acceleration (su- 
peracceleration) can only be provided in the context of 
EXQ through the time variation of Newton's constant 
G. A decreasing G{t) with time favors accelerating ex- 
pansion and allows for superacceleration. This physical 
argument is reflected in our results. We found that given 
Go — (neglecting gi since |5i| < 10"'' from solar system 
tests) a Go > ((72 > 0) decreases the forbidden sector 
and allows superacceleration {wo < —1). But Go ~ 
with Go > implies that we are currently close to a 
minimum of G{t) with G{t) being larger in the past ie 

^^l + lg,iHoit-to)f (3.1) 
Gro ^ 

Therefore a decreasing G{t) corresponds to Go > 
(92 > 0) (see Fig. 4) which in turn implies smaller forbid- 
den sectors and allows for superacceleration in agreement 
with the above physical argument. 

An additional interesting point is related to the con- 
struction of the contours of Figs. 1-3. The SNLS 
data analysis involved in the construction of these con- 
tours did not take into account the possible evolution of 
SnIa due to the evolving G(t). It is straightforward to 
take into account the evolution of G in the SnIa data 
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FIG. 4: A g2 > implies a decreasing Newton's constant 
~ 1 + \g2{Ho{t - to)f around t = to). 



analysis along the lines of Rcf. [38*1 . In order to test the 
sensitivity of these contours with respect to the evolu- 
tion of G we have repeated the contour construction 
assuming a varying G according to the ansatz 



G(z) =Go(l + a 



(1 



r) 



(3.2) 



which smoothly interpolates between the present value 
of G = Go and the high redshift value G = Go(l + a) 
implying that 

« = ^ (3.3) 
t^o 

The parametrization p.2p is consistent with both the so- 
lar system tests [s^l — 0) and the nucleosynthesis 
constraints 391 



1 G„ 



Go I 



Go 



< 0.2 



(3.4) 



at Icr, for |a| < 0.2. 

It is straightforward to evaluate gi, §2, <?3 hi terms of a 
using the parametrization (|3.2p and equation (|2.15p . The 
result is 



(3.5) 
(3.6) 
(3.7) 

(3.8) 

and repeated the SNLS data analysis taking into account 
the evolution of the Snia absolute magnitude M due to 
the evolving G [H, SO, Ei El as 



51 = 
92 = 2a 

33 = 3a (-1 + 3wo(-l + ^am)) 
We have considered the value of 
a = 0.2 



M = Mo + 



15, G 

t^^^g;; 



(3.9) 



The steps involved in this analysis may be summarized 
as follows: 

• Use the following magnitude redshift relation to fit 
to the SnIa data 



15 G(z) 

mth{z) = Mq + BlogdLiz) + —r^og—^ 

4 (jQ 



(3.10) 



where G(z) is given by p.2p and dL(z) is con- 
nected to H{z) in the usual geometrically defined 
way (|2:43l) ie 

dL{z)^{\ + z) j\z'^^ (3.11) 

• Use an H{z) parametrization that incorporates the 
evolution of G (given by (|3.2p ) in the context of the 
CPL parametrization (|2.32p ie 



H\z) 



G{z) 
Go 



Hl[no„^{l + zf + 



+ {l~ f^om)(l + z)3(i+"'°+«'i)eTiW](3.12) 
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FIG. 5: Taking into account a G variation consistent witii tlie 
nucleosyntiiesis and solar system bounds introduces modifica- 
tions to the Wo — wi contours which are not significant for 
the currently available Snia data (for the dashed line a — 0.2 
and the continuous a — 0). 



• Minimize the expression 



X {wo,wi) 



(3.13) 



and find the corresponding 
along the lines of Refs. [36 



la and 2a 
M 



contours 



The resulting contours in the wq —wi space are shown 
in Fig. 5 (dashed lines) superposed with the correspond- 
ing contours constructed by neglecting the evolution of 
G (continuous lines). The change of the best fit wq — wi 
values and of the corresponding errorbars is minor (espe- 
cially in the wq direction) and such that our main con- 
clusion regarding the limiting values of g2 remain practi- 
cally unchanged even after the evolution of G is taken 
into account in the analysis. This justifies neglecting 
the evolution of G in the construction of the con- 
tours for our purposes. However, as the quality of SnIa 



data improves, it becomes clear from Fig. 5 that the 
effects of an evolving G consistent with nucleosynthesis 
and solar system constraints on the data analysis can be 
significant! This possibility is further amplified if there 
are additional effects of an evolving G on the SnIa data 
analysis. For example [13,13] it is possible that the time 
scale stretch factor s involved in the SnIa data analysis 
[l], i36] and arising from opacity effects in the stellar at- 
mosphere may have a dependence on the Chandrasekhar 
mass and therefore on Newton's constant G. We have 
shown however that even if such effects are included in 
the SnIa data analysis our results of Fig. 5 (dashed line) 
do not change more than 10%. 

An interesting extension of the present work could in- 
volve the investigation of the limits of other dark energy 
models like braneworld models or barotropic fluid dark 
energy in the context of scalar-tensor theories (thus ex- 
tending the results of Scherrer [IHl)- The identification 
of the forbidden observational parameter regions for such 
models could be combined with the constraints of local 
gravitational experiments to test the internal consistency 
of these models. Alternatively the present work could 
be extended in the direction of finding the limits and 
cosmological consistency of specific classes of potentials 
in the context of EXQ thus also extending the work of 

Numerical Analysis: Our numerical analysis 
was performed using a Mathematica code available at 
|http: / /leandros. physics. uoi.gr/exqlim/exqlim. htm] 
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